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ABSTRACT 

The paradigm of the neutralino dark matter predicts that the first gravitationally bound objects are 
earth-mass sized microhaloes, which would emit annihilation gamma-rays. Here we show that, though 
the flux from individual nearest microhaloes is extremely difficult to detect, meaningful constraints 
on their survival probability and internal density profile can be set by requiring that the galactic and 
extragalactic gamma-ray background flux from the microhaloes does not exceed the existing EGRET 
background data. Possible disruption of microhaloes by stellar encounters does not significantly reduce 
the background flux. If the probability for microhaloes to survive the hierarchical clustering process 
of dark matter is as large as indicated by a recent simulation, they could be a significant component 
of the observed background flux in some photon energy range, even with the standard annihilation 
cross section and conservative internal density profile of microhaloes. The integrated gamma-ray flux 
from microhaloes in the halo of the Andromeda galaxy may also be detectable by observations in the 
near future. 

Subject headings: dark matter — gamma rays: theory — Galaxy: halo — galaxies: haloes 



1. INTRODUCTION 

The neutralino predicted by the supersymmetry the- 
ory of particle physics is the most promising candidate 
of the cold dark matter (CDM). The theory predicts 
that neutralinos should annihilate and produce high en- 
ergy particles such as gamma-rays, and the detectabil- 
ity of these particles from dense regions like the Galac- 
tic center (GC) has been discussed intensively in the 
literature. [See iBertone et alJ i|2004l) for a review and 
references therein for the earlier papers.] The struc- 
ture formation theory predicts hierarchical substructure 
in a dark matter halo, and this dumpiness is expected 
to enhance the annihilation signal. Previous investi- 
gations, however, considered substructures only to the 
first order (i.e., isolated subhaloes in a halo) and masses 
larger than ~ 10 6 M Q , mainly because of the limitation 
of cosmological N-body simulations lUllio et alJ l2002t 
ITavlor and Silk 20MM 

sacss er and M annheim 2005|). 
Recently attention has been turned to substructures 
on much smaller scale, especially the first gravitation- 
ally bound objects in the cosmological evolution. If 
the dark matter is the neutralino, any density fluctu- 
ations of mass scales smaller than M m h ~ 1O -6 M0 
are washed out by collisional damping and subsequent 
free streaming in the early universe. Then the first 
objects form as "microhaloes" with a mass ~ M m h at 
z ~ z n \, where z n \ = 60 ± 20 is the epoch when the rms 
linear density f luctuation at this mass scale, a(M m h), 
becomes unity lIHofmann et al.N2001l IGreen et ;iljl2()()i 



l2tfjilBerezinskv et al.l2003HLoeb fc Zaldarriaea Il2005jl 
Therefore a considerable part of the neutralino dark mat- 
ter should have collapsed into these earth-mass objects, 
at least once in the cosmic history. 
However, it is highly uncertain how much fraction of 
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these microhaloes can survive the subs equent hierarchi- 
cal str ucture formation until present. IBerezinskv et all 
(2003) estimated that only 0.1-0.5 % of microhaloes sur- 
vive, because of tidal disruption when they are taken 
into larger haloes. This is, however, a completely an- 
alytic estim ate and uncerta i nty m ust be large. On the 
other hand, iDiemand et al.l l)2005|) . based on a N-body 
simulation, argued that about 50% of the total halo 
mass Mtot is in the form of substructure with a sub- 
halo mass function dN/dM oc in the mass range 



io- 



< M < 1O 1U M , where \i 



This indi- 



cates that the mass fraction of microhaloes is at least 
Ml h [dN(M mh )/dM]/M tot ~ 1.3%. Furthermore, the 
nested nature of hierarchical structure formation pre- 
dicts that microhaloes may be embedded in larger sub- 
haloes, which may again be embedded in even larger 
ones. Counting microhaloes in larger mass subhaloes 
up to M ~ 10 10 Mq will further increase the tru e num- 
ber of earth- mass microhaloes. IDiemand et alJ (2005) 
used 500 pc -3 as the number density of such micro- 
haloes in the solar neighborhood, which is about 7% of 
the standard dark matte r density at the Sun' s location, 
p(Ro) = 0.3 GeV cm" 3 ijBertone et al.l l20041. In addi- 
tion to the tidal disruption by hierarchical clustering, mi- 
crohaloes may also be dest royed by tidal interaction with 
stars in the Galactic disk (jZhao et alJl2005l iMoore et alJ 
2005), but again estimates are controversial. 

It is obvious that such microhaloes could have sig- 
nificant impact on the detectability of the annihilation 
signal. Here we show that both the galactic and ex- 
tragalactic gamma-ray background radiation (hereafter 
GGRB and EGRB, respectively) give a strong constraint 
on the existence of such microhaloes. Since the micro- 
halo survival probability is highly uncertain, we simply 
parametrize this quantity as / SU rv, and try to estimate 
how much constraints can be set from existing and future 
observations. It has been argued that neutralino annihi- 
lation cannot be a significant component of the observed 
EGRB data since it would overpredict the gamma-ray 
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flux from the GC beyond the observational upper bound 
ijAndd 12005ft . This argument does not apply here, be- 
cause the microhaloes are likely disrupted by strong tidal 
forces of the GC gravity fiel d and/or interaction with 
stars within ~ kpc of the GC l)Diemand et al]l2005ft and 
hence the visibility of the GC is not enhanced by the 
microhaloes. However, the mass included within 1 (10) 
kpc i s only 0.3 (6) % of the total mass of the Galactic 
halo !|Klvnin et all 12002ft . Therefore the total gamma- 
ray flux from microhaloes in a galactic halo, to which 
the EGRB is related, is hardly affected by tidal disrup- 
tion in its central region, if microhaloes trace the mass 
distribution. 

2. GAMMA- RAY BACKGROUND FROM MICROHALOES 
2.1. Estimating GGRB and EGRB Flux 

In this work we assume that microhaloes are formed 
at redshift ~ z n \, and a fraction (1 — / surv ) of them are 
immediately destroyed by tidal forces in subsequent hier- 
archical structure formation. After that, we assume that 
the number of microhaloes and their density profile are 
kept constant until present, except in regions very close 
to galactic halo centers. These are not unreasonable, 
since the tidal force within an object is proportional to its 
internal density, and the mean density within virialized 
objects decreases with the cosmic expansion in propor- 
tion to the mean background density. Therefore we ex- 
pect that the tidal disruption at a typical location within 
a dark halo should occur most efficiently at redshift not 
very different from z n \, when isolated microhaloes are 
first taken into larger objects. 

Consider a region in the universe with the mass 
scale M m h whose linear fractional overdensity is 5 = 
5p/p oc (1 + z)^ 1 . According to the standard struc- 
ture formation theory, the height over the rms, v = 
5/cr(M m h), obeys to the Gaussian, i.e., the comoving 
number density of microhaloes given by dn m \yjdv — 

./snrvfivPm-t..n ex p(-i/ 2 /2)/(y / 27rMrr| h), where we use the 
WMAP values i|Snergel etafl 12003ft for the present-day 
critical density p C rit,o and the neutralino density param- 
eter Q x — 0.22. At z ~ z n \ the universe is flat and 
matter-dominated, and the region will collapse and viri- 
alize when 6 grows to S c = 1.686 at z = z v j r , where 
(1 + 2 v i r ) = (1 + z u \)v/8 c . The internal density of mi- 
crohaloes is given as p e &(v) — 187r 2 / c (l + z V ir) 3 ^x/>crit,o- 
Here, f c is the enhancement factor from the virial density, 
to take into account the density profile of each microhalo; 
we found that the mass-weighted mean density, which 
is proportional to the annihilation rate, is increased by 
f c = 6.2 for the microhalo profile 3 found in the simula- 
tion Ipiemand et al.ll2005ft . 

Then, integrating over is, the comoving annihilation 
rate density is given by: 



J<f v 



M mh p e fi(is) 



2m 2 



dn r 



dv 



dv 



(1) 



where m x is the neutralino mass and (o~ xx v) being the 
mean velocity-multiplied cross section of neutralino an- 
nihilation. Throughout this letter we use the standar d 



value of (o\ 



3x 10~ 26 cnrV 1 l|Bertone et all2004ft . 



3 The a/37-profile with (a, /3, 7) 
tration parameter c = 1.6. 



(1, 3, 1.2) and the concen- 



and scaling for different values is o bvious. The possibl e 
range of m x is ~ 30 GeV - 10 TeV ijBertone et al.ll200^1 . 
and we use m x = 100 GeV for calculations below, un- 
less otherwise stated. Since p c s oc ^ 3 , the annihilation 
signal comes mainly from microhaloes of ~ v p sigma 
fluctuation, where v p = V3- Microhaloes with v <C 1 
may not gravitationally collapse or would be disrupted 
by subsequent structure formation, but the integration 
is not very sensitive to the lower bound, and hence we 
take vi — 0. Note that, by this formulation, / SU rv is ef- 
fectively the survival probability for microhaloes having 
relatively high density fluctuation of v ~ v p , and their 
mass fraction in the total dark matter, f m , is related 
as: f m ~ (l/\/27r)exp(-f 2 /2)/ surv ~ 0.09/ surv . Assum- 
ing that the spatial distribution of microhaloes traces 
the smoothed mass over larger scales, the number den- 
sity of microhaloes with v <~ v v around the solar sys- 
tem is n{R ) ~ f m p{R&)/M mh ~ 680/ surv pc" 3 . This 
number is cl ose to the estimate ba sed on the simulation, 
~ 500 pc -3 ijDiemand et alJl2005^ . indicating that / surv 
could be of order unity. 

Now we can calculate the EGRB photon flux from mi- 
crohaloes per steradian, as: 



dF 7 
dEL 



4tt 



r» l ^dt dn 7 [{l + z)E 1 ] 
J dz dE-y 



+ z) , (2) 



where E 7 is the gamma-ray energy and its spectrum pro- 
duced by an annihilation, dn 7 /dE 7 , is calculated using 
the an alytical fitting formula given in iBergstom et alJ 
l)2001tK (We consider only the continuum gamma-rays.) 
The integration up to the redshift z n \ is an approxima- 
tion, but it is almost insensitive to this upper bound 
since annihilation at z < 1 is dominant to the EGRB. 
High energy gamma-rays may be absorbed by interac- 
tion with the cosmic infrared background. However, in 
this letter we consider only E 1 < 100 GeV, and at this 
photon ener gy the optical depth become s unity only be- 
yond z ~ 2 ijTotani and Takeuchll2002jl . Therefore our 
result is hardly affected by the absorption. 

Next we calculate the GGRB flux. The annihilation 
rate per unit dark matter mass in a region smoothed 
over larger scales than microhaloes is given by N xx — 

A/" xx /(O x p cr it,o)- Then we obtain the GGRB flux per 
steradian as: 



dF~ ( N xx dn 1 {E 1 ) 



dE~, 



4tt dE. 



1 



dl 



(3) 



where the integration is over the line of sight, for the 
smoothed dark matter density in the Galactic halo, p sm . 
It should be noted that the flux is proportional to the 
line-of-sight integration of Pg m , not p^ m as in the case of 
diffuse matter distribu tion. We use two sp herically sym- 
metric models of p sm lIKlvpin e t al.l 120 02ft: one ha s the 
Navarro- Frenk- White fNFW. iNavarro et alJ 11996^1 pro- 
file but the other is modified from the NFW profile by 
adiabatic compression of dark matter responding to the 
baryon infall. As discussed above, microhaloes are ex- 
pected to be destroyed by tidal forces in the inner region 
around the GC. Therefore, as a simple model, we intro- 
duce the disruption radius Rd within which microhaloes 
are completely disrupted, while they are all preserved 
outside Rd with an abundance proportional to / surv . 
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2.2. Comparison to Observations 

The calculated EGRB flux is shown in Fig. 2] for m x = 
100 GeV and 1 TeV. If / surv ^0.1, the microhaloes make 
a significant contrib ution to the observed EGRB flux 
{Strong et alJl2004a|) in some photon energy range. The 
pre dicted flux is much higher than those in earlier stud - 
ies (|Ullio et alJl200l iTavlor and Silkll200l IAnddl2005j) : 
this is mainly because the microhaloes formed much 
earlier than galactic haloes considered in these studies 
(M > 1O 5_6 M0) and hence have much higher internal 
density. The GGRB flux from the microhaloes, which is 
the mean in the all sky except for the disk region 4 , is 
also shown in this figure. It is found to be comparable to 
the EGRB, indicating that it could also be a significant 
component of the observed flux. It should be noted that, 
since the "observed" EGRB data was estimated as the 
residual after the subtraction of the cosmic-ray interac- 
tion model in our Galaxy, not only the EGRB but also 
the GGRB expected from microhaloes should be com- 
pared to the observed EGRB data. 

Figure [3 shows the predicted total flux (EGRB + 
GGRB) as a function of the Galactic longitude. It can 
be seen that, if Rd > 5 kpc, the anisotropy of the back- 
ground flux is at most a factor of 2. This is acceptable, 
considering the precision of the EGRB flux measure- 
ments and possible systematic uncertainties in the fore- 
ground subtraction llSreekumar et al]ll998l iStrong et al.l 
l2004aHKeshet et al.ll2004ft . Even smaller R d may also be 
allowed, since the anisotropy close to the GC would be 
hidden by strong background flux from cosmic-ray inter- 
actions in the Galactic disk. 

In fact, evidence for a diffuse gamma-ray halo towards 
the GC that cannot be explained by the stan dard cosmic- 
ray i nteraction model has been reported {Dixon et alJ 
1998), at a flux level similar to the EGRB. This gamma- 
ray halo might be explained by the microhaloes with an 
appropriate choice of Rd- The GGRB sky distribution 
is expected to show strong small-scale anisotropy by the 
complicated substructures in the Galactic halo. Since 
annihilation signal peaks in rather narrow photon energy 
range, the gamma-ray energy dependence of the GGRB 
anisotropy may be used to examine the contribution from 
microhaloes. 

The EGRET background data around the disk region 
shows so-called GeV excess over the standard prediction 
from cosmic-ray interaction, which is about 1-2 orders 
of m agnitude higher than the EGRB along the Galactic 
disk {Hunter et al.lll997j) . The flux level of the excess 
might be achieved by microhaloes with a large boost fac- 
tor from our prediction above, which is, in fact, not very 
unlikely (see below). However, still the GeV excess seems 
difficult to explain by the microhaloes, because the spa- 
tial distribution of the GeV excess is clearly associated 
with the Galactic disk while distribution of microhaloes 
is expected to be more spherical. Note that the GeV ex- 
cess can also be explai ned by modification of the cosmic- 
ray in teraction models {Strong et al.l2004bllKamae et alJ 
2005). 

It should be noted that the internal density profile of 
the microhaloes used above is conservative in a sense that 

4 Th e removed disk region is the same as defined bv lStrong et al.| 
12004a), for a consistent comparison between the prediction and 
the data. 



it predicts relatively low annihilation luminosity. Though 
we have used the aff'y profile w ith 7 = 1.2 following the 
fitting bv lDiemand et all {2005) (p oc r'^ with r -> 0), 
their simulated microhaloes have 7 ~ 1.7 to the resolu- 
tion limit of the simulation. They also noticed the sim- 
ilarity between their simulated microhaloes and galac- 
tic haloes shortly after the formation or major mergers, 
showing a single power law profile with slopes of 7 ~ 1.5- 
2. If 7 > 1.5, annihilation luminosity diverges in the cen- 
ter, and assuming that the maximum density is limited 
by annihilation time scale, p max (a xx v) / (2m x ) < (10 10 
yr)" 1 , we find that the enhancement factor f c = 6.2 used 
above (for 7 = 1.2) will be boosted up to f c — 31, 190, 
and 1.4 x 10 4 for 7 = 1.5, 1.7, and 2.0, respectively. Here 
we used consistently (and conservatively) a low value for 
the concen tration parameter, c = 1.6, as found in the 
simulation {Diemand et alJfeOOSlK 

Provided that annihilation cross section is close to the 
standard value, the two major astrophysically uncertain 
parameters are / surv and f c . We note that the back- 
ground flux is only weakly dependent on M m h. The 
formation redshift z n i depends on M m h, but only very 
weakly, because c(M) is only weakly dependent on M 
in small scales under the standard CDM density fluctu- 
ation spectrum and hence fluctuations over a wide range 
of mass scales will become non-linear at similar redshifts. 
Hence we calculate the excluded region in the / SU rv-/c 
plane in Fig. [21 by requiring that the predicted flux does 
not exceed any of the observed EGRB data. 

3. DISCUSSION AND CONCLUSIONS 

The observability of the nearest microhalo from the 
Earth is of great interest. The expected photon flux is 
given by: 

F Y 7 M mhPcS (o- xx v) , . 

~ 4tt [n(i? )]" 2 / 3 2m 2 , 1 > 

= 1.5 x 10- u /su / ^(/c/6.2)y 4 om^ 2 cm^s" 1 , (5) 

where rri2 — m x J (100 GeV) and V 7 = 4OY40 is the 
photon number yield from one annihilation; typically 
one annihilation produces 30-50 continuum gamma- 
rays and more than 80% of them are above 100 MeV 
{Gondolo et alJl2004t) . Unfortunately this flux is much 
smaller than the point source sensitivity of the EGRET, 
and even of the future GLAST mission, and hence detec- 
tion of individual microhaloes is unlikely. It is impossi- 
ble to greatly enhance the detectability by boosting up 
the internal density factor f c , because it would seriously 
overpredict the background flux far beyond the observed 
level. 

The detectability of gamma-rays from microhaloes 
in the nearby extragalactic objects is also intrigu- 
ing. Here we estimate the flux expected from M31, 
d = 770 kpc from the Earth. The dark matter 
mass enclosed within 13 kpc, corresponding to the 
position accurac y of EGRET (~ 1°), is M DM ~ 
1.5 x 10 n M o {KlvDin et al. 2002). Then we ex- 
pect photon flux of Y 1 N xx M DM /{4.ircP) ~ 7.7 x 
10~ 9 /surv(/ c /6.2)l4om^ 2 cm~ 2 s _1 , which is interest- 
ingly very close t o the EGRET upper bound (1.6 x 10~ 8 , 
IBlom et al.lll999h . indicating that there is a good chance 
of detecting annihilation flux from microhaloes in M31 
by the GLAST mission or next generation air Cerenkov 
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E r [GeV] 

Fig. 1. — The background gamma-ray flux from neutralino annihilation in the microhaloes. The Galactic (dashed), extragalactic (dot- 
dashed), and the total (solid) components are shown. The two cases of m x = 100 GeV and 1 TeV are presented, with / surv = 0.35 and 
1, respectively. The internal density profile parameter f c = 6.2 is conservatively assumed. The baryon compressed NFW profile for the 
Galactic halo and R d = 5 kpc are used for the Galactic component. 



telescopes (ACTs), with a larger m x favoring the latter. 
In contrast to the flux from the center of cuspy density 
profile, we expect rather diffuse flux distribution around 
the center of M3I, because microhaloes that are very 
close to the M31 center are likely disrupted. The im- 
proved angular resolution of the GLAST or ACTs will 
be able to resolve it. 

In conclusion, annihilation gamma-rays from the mi- 
crohaloes could be a significant component of the ob- 
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Fig. 2. — The longitudinal distribution of the background flux (GGRB+EGRB), assuming m x = 100 GeV, / sur v = 1, and f c = 6.2. 
Three different values of disruption radius R^ are used as indicated. Two different density profiles of the Galactic halo are used: the 
baryon-compressed NFW (upper thick curves) and the original NFW (lower thin curves). The level of the predicted isotropic EGRB from 
microhaloes is also indicated. 




fc 

Fig. 3. — The excluded region in the space of the two parameters, the survival probability (/ surv ) and the enhancement factor by 
the internal density profile of microhaloes (/ c ). The microhalo mass fraction f m in the total dark matter mass is related to /surv as 
fm ~ 0.09/ S urv Several curves are depicted for different neutralino masses (m x ) as indicated, and the upper-right regions are excluded 
because the predicted background flux will exceed the observed data. The compressed NFW profile for the Milky Way halo and Rj = 5 
kpc are assumed for the GGRB component. The values of f c corresponding to several values of the inner slope index (7) of internal density 
profile of microhaloes are marked by vertical thin solid lines. 



